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The process of combining telemetry signals received at multiple antennas com- 
monly referred to as arraying, can be used to improve communication link perfor- 

muhtelTr ^ f PaCe NetWOrk < DSN )■ B y coherently adding telemetry from 

?SNR ove th7 » signal- to- noise ratio 

Lt f h h r achwvable mth Hny sin S lc antenna in the array. A number ofdiffer- 

pJliteraZl STtT , been , pr ° posed and thep Performances analyzed in 
past Merature [1,2]. These analyses have compared different arraying schemes un- 

and th TT 1 t m tbC , tilgnals COntain additive white Gaussian noise (AWGN) 
and that the noise observed at distinct antennas is independent. 

in Whele ™ ^wanted background body is visible to multiple antennas 

m the array however, the assumption of independent noises is no longer applicable 
A planet with significant radiation emissions in the frequency band of interest can 

touTof SOU ^ c of ^related noise. For example, during much of Galileo’s 
upiter, the planet will contribute significantly to the total system noise at 
ions ground stations. This article analyzes the effects of correlated noise on two 

^Sf^ScL C T nt7y , being Ti Mered f ° r DSN appIicati ons: full-spectrum 
combining (FSC) and complex-symbol combining (CSC). A framework is presented 

for characterizing the correlated noise based on physical parameters, and the impact 

of the noise correlation on the array performance is assessed for each scheme. 


I. Introduction 

may be used to eohetently track signals that are too weak to te”r^ by“ an^^r 

processes that are used to combine and Site ' he signal ^nH t" f Sy ” Chr0 “ k “ 1 °" 

" hich is a 
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Previous analyses E 

used an additive white Gaussian noise (AWGN) model to However, if a strong radio 

assumed the noise waveforms received at ^tmct antenn^ d P observations at different 

=r^:”/cr s ^ 

^ m P «“. h Srd “foHsTr^Xr'ralysls is thus needed to characterise the performance of 
arraying schemes in cases where correlated noise is present. 

Prior work has been conducted on this subject but has phy^ai 

by Dewey [3] examines correlated noise effects 1 * ° £ lnto a ,,. ol][lt properties of the source and 
considerations. A correlated noise model is pre ’ arrayed symbol SNR relative to a case 

the array geometry The impact of the - observation of the Galileo 

of uncorrelated noise ls.then analyz . Australia complex However, Dewey’s study does not 

spacecraft from a synchronization in telemetry arraying, which are dependent on 

used. Thus, the analysis does not identify the relative advantages and 
disadvantages of different arraying schemes under conditions of correlated nois . 

The purpose of this article is t0 

(FSC) and complex-symbol combining (CSC) array g svstems is presented. Parameters used 

to understand the physics underlying background I "o.se >n I ^“^nd explained. Sections III and IV then 

“ g h :,Tand S “re n pIrle,s £ nZd t ^uLT^rformance of both arraying schemes in this scenarro. 
Finally, Section VI summarizes the mam results of the work. 


II. Background Noise Properties 


ttaCKgrounu nviac * ■ * — 

Here we present basic 

remainder of the analysis. 1 he discussion treatment of the subject material may refer to a 

performed by Dewey alluded to earlier ,3,. 

Consider hrst the effect o, a b^roun, 

an antenna consists of both therma noise tvnirallv have an emission spectrum that varies 

radio sources in the antenna's field o, The Wndwidth of interest’ The 

brightness distribution over the antenna’s reception pattern, i.e., 


T s = ^ J J B(s)Pn(s) ds 


( 1 ) 


where A e is the effective receiving area of the /hJ/sf (tr^tands for steradian, a measure of solid 

S® — — — “ — “ 

2 Ibid. 
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The one-sided power spectral density of the noise due to the source is then given by N - kT NntP ^ 

zziz: r ce ls c — 1 in the peak » f - ^ 


t s = 



(2) 



( 3 ) 


where S is the total flux density of the source in W im^ I H ? Ac i , . i 

and the spacecraft increases, the background source loves on, tf 

~L“ 


^n„n t (r) = E[h t (t)n* k (t - T )] = a (4) 

:«S 


R n^h k (r) = a S(t) 

It can be shown [3,4] that the cross-power spectral density level is given by 

/ J B(i)y/P Ni (S)P Nk (i) e i»'f'>B»-*/c dS 


( 5 ) 


~^ek 

a — — *- — - 

2 




-\V\e 


j<t>v 


( 6 ) 


where f Q is the observation frequency, and c is the sneed of liaht q v m8 / T . 

applications, the gnuntity |K|e«- is kuowu as the complex Liility £ th/LL A ^2 

a^i daT^I 8 E<1 ' i 6> “ e made ^ Fi ' 5t ' ” K th “ th * ex P° ne ntial term e>-/-ft. “ 

y //o 2 /c r> , ere -B,* p is the projected baseline length in the direction of the smirrp Tf u 

being the angular radius of the source) Irvl — ► n a nH fi, • u . ‘ lkj> ^ c /(/o#s), with 

uncorrelated. By contrast for B ”T ’iH p non* observations due to the source become 

achieves its uppeHimtTamdy " ^ ^ the spiral density 
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( 7 ) 


a 


\J A €t A e 


the greater the impact of a background body on the array. 

Finally, we introduce the correlation coefficient, describing the degree of correlation that exists between 
the noise at two antennas, defined as 


A 

Pik = 


\/ N 0t N 0k 


( 8 ) 


Note that in the upper limit (i.e., source size 
becomes 


small compared to fringe period), the correlation coefficient 

( 9 ) 


Pik 


ft jga 
V T t T k 


n T T are the source temperatures at antennas i, fc, and T), T k are the total system temperatures 
ll “Si ^ Thus, the greater the contribution of the source to the total system temperature, 

the higher the correlation coefficient, as is intuitively expected. 

Combining Erp. (5), (6), and (8), the cross-correlation function for the noise observed at two antennas 
can be expressed as 

Ra„n k {r) = /WiVo.WoT k 6 ( t ) (10) 

where is used to express the correlation phase, denoted by <£„ in Eq. (6). 


Ill Full-Spectrum Combining Performance 

Given a mathematics, description of noise ^ “mmS 

here briefly. Assume the ^ f downconverted 

antenna (i.e., the antenna with the high / •) . , • w hi c h can be thought of as a 

to baseband 3 by local oscillators in phase qua ra ure ' compensate for differing arrival times 

single complex signal, is then shifted in time y s °™ e * baseband signals are then aligned in phase, 

the si8 " al 15 processed by 8 

single carrier, subcarrier, and symbol loop. 

Two quantities used to describe .arraying tsW ratto of the ideal 

“mbffiSfol 2? J55?S- to the id/., symbol SNR of — 1 W- Here, H-P meaus the 


‘ cult, we 

not dependent on what frequency processing is done at. 
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Fig. 1. Full-spectrum combining. 


lt (u - per “ 

of which arraying technique is used, since synchronization losses are ig"3” 8 £s O “d Ttf 

ssssssss 

f ,l n] f , ,f atl ° n for tull - s Pectrum combining and complex-symbol combining was commited in 

A. Ideal Arraying Gain 

2/iW =Si(0+«i(0 


- v/2-Pr, cos (w [F t + 6 t + Ad(t) 




\/2Pc, COS (w/F< + 0*) y/2P Di d(t) sqr (u sc t + 6 SC ) sin(w /F t 4- 9 t ) + m(t) (11 

in radians/s; 0 SC is the subcarrier phase in radians- P i thf ’ " h subcarrier frequency 

Pt cos 2 A- P,a ic ihnrii radians ’ P c , is the carrier power in watts, given bv Pr = 

T. cos A, P D , » the data power watts, given by P D , , P% sin 2 A; and „,«) is an addije white 
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Gaussian notes process with one-sided power spectra, density W/H, The corresponding complex 
baseband signal is given by 

( 12 ) 


y,(f) =- Si(t) + hi(t) 

= s/Pc i e i{ ' Jbt+ei) + jy/PoAt) sqr(w 3C t + 6> sc )e' (Wbt+< ' ,) + n<(t) 


(13) 


2SL3 no£ 

The spectrum of the baseband telemetry is shown in Fig. 2. 

f b = baseband carrier frequency 

f sc = SUBCARRIER FREQUENCY 



Note that the bandwidth needed 

determined by the subcarrier frequen y, fsc ad ’ , * transmits portions of the 

AS an alternative ,0 the method Each signal can 

modulation; the total tr^sm^ssnm^^tandwidt^s^ J ^ ^ „ ttat the processing required 

is mentioned briefly in (l). no% e , , . l modified for each mission. For 

SmpC- “if' ^^s b o7t"bSc implementation of FSC de^ribed in H . ■-«•« * — « h “ 
a more bandwidth-economizing option also exists. 

Let the phase different tetween the phmmomm^by^an^tirnate ofthis ^^tity) 0 io tl^idi^^lwnt 
Tilted s[gi are ,he P „ multiplied by prespecihed weighting factors, ft, and summed. 
The combined signal is thus given by 


Vcombit ) = S co mb{t) + «comb(f) 


(14) 


1=1 1=1 


(15) 


_ Y2p %e >* u +J\/PD. d ( t ) sqr {u) ac t + O sc 


^ e j(te'b t + 0 ‘)^ 


i=l 


+ 2A^ ,lfl i ( t ) 


(16) 


1=1 


216 


where the weights ft are chosen to satisfy the condition 


ft = 


Fh Noi 

Pt x N oi 


(17) 


for i = 1, ■ • • , L. It is shown in [1] that these weights maximize the combined SNR when the noises 
ni(£) are independent. Note that this is not necessarily the optimal choice of weights for the correlated 
noise case, as pointed out in [3], Furthermore, the optimal choice of phases used to array the signals is 
not necessarily the relative signal phases, <pu. Using the phases <pu will certainly maximize the arrayed 
signal power, but not necessarily the ratio of signal to noise power, which is the relevant criteria for 
optimization. The problem of optimal combining weights and phases for signals with correlated noise 
has been analyzed in [5], where the results are applied to an array of antenna feed elements. However, 
computation of these weights requires knowledge of the pairwise correlations between the noises, aye^S, 
for all i,j pairs. A scheme can be devised to estimate the required parameters in real time and modify the 
weights accordingly, but would significantly complicate the problem. Our goal, instead, is to determine 
the performance impact of the correlated noise assuming the traditional combining scheme is used. 

The total combined signal power, Pt, is given by 

Pt = E [5 com &(t)] E [s* om b(£)] (18) 

If the relative signal phases are estimated perfectly (i.e., <pu = <p\ z for i = 2, • ■ • , L), the combined signal 
power becomes 


P T = P Tl 


( L L L 

I i - 1 1=1 i=i 

\ 


(19) 


where 7l = [(PrJAPr, )][(N 0l )/(N 0i )]. 

The one-sided power spectral density of the real and imaginary parts of the combined noise is given 
by 


~ 2 ^ Comb ^ ™comb(t)] ( 20 ) 

where B is the one-sided bandwidth of the noise waveforms. Note that the factor of two in the denominator 
of Eq. (20) results from the fact that the real and imaginary parts of the noise each has half the power 
of the complex noise. From the definitions of power spectral density and cross-power spectral density, it 
follows that 


E[h t (t)n*(t)} 

= 2N 0 ,B 

(21) 

E [n t (i)n«(i)] 

= 2 a ije j ^B 

(22) 


Equations (20), (21), and (22) can be combined to find the power spectral density of the combined noise, 
yielding 
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(23) 


N 0 = N ( 


0 1 


L L 


e ^ + n s 4i ' i) 


»= 1 i= 1 i=! 


The P T /JV 0 of the combined signal is thus given by 


Pt _ ft-, (^i 7i ) 

N 01 ^f =1 7 i + Ef=i Ef=j \flPTi Pij 


(24) 


where ip,j = 4>? ~4>ij- The parameters p tJ and ip i3 describe the relevant statistics for the noise correlations 
between the various antenna pairs and determine the correlated noise impact on the ideal arraying gain. 

The combined signal is finally processed by a single carrier, subcarrier, and symbol loop. Assuming 
perfect references at each of these three stages, the symbol SNR of the arrayed signal becomes 

(Z) L -i7*) 2 

Ef=i 7, + £f=i *** 


Ga (25) 


SN Rideal 


2P Dl 

N oi Rsym 


2Pp l 
N Cl Rsym 


where G a is the ideal arraying gain due to combining the signals. Note that setting all the noise correlation 
coefficients p i3 to zero results in G A = E,= i 7*. which is the ideal arraying gain in the case of uncorrelated 
noises, as discussed in [1]. 

Further note that the ideal arraying gain in the presence of correlated noise can be higher or lower than 
the uncorrelated noise case, depending on the phases ipi j- This point can be understood by considering 
an array of two equal antennas (i.e., 71 = 72 = 1.) Figure 3 shows values for G A for two equal antennas 
as a function of p and ip. For p = 0, the ideal arraying gain is a constant 3 dB, as expected. Now suppose 
the noises have some nonzero correlation coefficient, p, and some correlation phase, <p n . If ip = 0 deg, 
then the phase difference of the spacecraft signal as observed by antennas 1 and 2, <p, is equal to the noise 
correlation phase <p n . Thus, phase aligning the two signals also phase aligns the correlated component 
of the noise. The noise from the background source adds maximally in phase, and the combined noise 
power increases. Thus, the combined SNR decreases, and hence the arraying gain falls below 3 dB. By 
contrast, if ip = 180 deg, phase aligning the signal results in combining the correlated component of the 
noise 180 deg out of phase. Thus, the noise combines destructively in this case, and the arraying gain 
is now greater than 3 dB. For intermediate values of ip, the arraying gain varies continuously from its 
minimum value at ip = 0 deg to its maximum at ip = 180 deg. 

B. Symbol SNR Degradation 

In practice, perfect phase alignment and ideal carrier, subcarrier, and symbol references are not avail- 
able. Some degradation in the arrayed symbol SNR is, therefore, incurred due to synchronization errors. 
To quantify the degradation, we first find the set of density functions for the phase alignment errors 
Atpu = 0n - (pu, i = 2, ■ • • , L. This set of functions is then used to compute the Pt/No of the arrayed 
signal. Adding in losses due to carrier, subcarrier, and symbol tracking, the symbol SNR at the matched- 
filter output can be computed. Finally, comparing the actual symbol SNR to the ideal symbol SNR given 
by Eq. (25) yields the degradation for full-spectrum combining. 
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1. Antenna Phasing. A set of phase estimates <j>u for i = 2, • ■ • , L are needed to align signals 
2 j'"\ L with signal L In the description of FSC given in [2], the phase difference between s x {t) and 
si(t) is estimated by filtering the two signals to some lowpass bandwidth B lp Hz, multiplying them, and 
averaging their product over T corr s. The phase of this complex quantity is then computed by taking the 

inverse tangent of the ratio of the imaginary to real parts. A block diagram of this scheme is shown in 
Fig. 4. 

The complex product of the baseband signals after averaging, Z, is given by 


Z 



+ fHpi W)(^7 Pl (0 + ^(0) * 


= WPc.Pc, + \/P Dl PD t H)e j *" 


P rp I (^5,n(0 4- nip x (t)fll p t (f)) dt 
corr J 


(26) 


where H is given by 


H = 



1 

2 P 


(27) 


and M is the highest harmonic of the subcarrier passed by the lowpass filter. The term h sn (t) is 
composed of signal-noise terms in the product and has zero mean. Note, however, that the noise-noise 
term, h/ Pi (t)h iPt (f )* , does not necessarily have zero mean, due to a possible correlation that exists between 
the two noise waveforms. The expected value of this noise product can easily be computed from the cross- 
power spectral density of fii(f) and hi(t); thus, 

E \ z \ = (V E c, Pc, + s/PdJPd^H)^ 1 * + 2pi t sj N 0l N 0i B ip e (28) 

Since <f>\\ is not necessarily equal to 4> u , the noise product introduces a “bias” to the estimate of the 
relative signal phase. This situation is represented pictorially in Fig. 5. The complex quantity E[z] can 
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rj(t) 


4(0 



Fig. 4. Conventional phase estimator. 



be thought of as a vector sum of a signal-to-signal correlation, 
Note how the presence of the noise vector biases the measurement 
The relative magnitude of these vectors is given by 


S, and a noise- to-noise correlation, N . 
of the phase of the complex correlation. 


A 

|iV| 



(29) 


For typical parameters, even relatively modest levels of noise correlation can lead to a substantial biasing 
effect in estimating the relative signal phase. For example, consider correlating two signals, each having 
& Pt/N of 20 dB-Hz with a 1-kHz correlation bandwidth. Even if all subcarrier harmonics are induded 
in the correlation, making H = 1, a correlation coefficient as low as p = 0.1 makes the ratiom Eq. (29) 
equal to 0 5. The phase estimates are then influenced more by the relative noise phases j6 lt than t e 
desired quantities <Pn, leading to a high amount of degradation in combining the signals. If the alternative 
method described in Section III. A is used, where the subcarrier harmonics are filtered individually prior to 
combining, the effective correlation bandwidth can be lessened, thus reducing the impact of the noise bias. 
Nevertheless, a practical implementation of full-spectrum combining requires a modified phase estimation 
algorithm if correlation levels encountered will generate significant biases. 


The method of phase estimation shown in Fig. 6 can be used for this purpose. Here each signal is 
filtered to some bandpass bandwidth B bp , and an additional complex correlation is performed between 
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Fig. 6. Modified phase estimator. 

the resulting waveforms. The center frequency of this filter is chosen so as to not capture any energy 
from the telemetry; this can be accomplished by locating the filter at an even multiple of the subcarrier 
frequency for example. After scaling the noise-only correlation by the ratio of the lowpass-tobandpass 
bandwjdths, this quantity provides an estimate of the contribution of the noise to the total correlation 
The bandpass correlation ^can then be subtracted from the lowpass correlation to compensate for the 
mean correlation vector |TV|. The compensated correlation can thus be expressed as 

Z = W Pc ' Pc - + y/PDi p D t H) + -L- f (n, iB (t) + h lpi (t)ht (t)) dt 

/ J-corr J Px ' 

B 1 f 

d r r > I ( t ) dt 

J ->bp corr J 


= (V P C\Pc\ + y/PDiPD.H)**" + N 


where the the noise term N now has zero mean. The phase estimate is then found by taking the inverse 
tangent of the ratio of the imaginary-to-real part of Eq. (30), i.e., 


4>u = tan 1 


W p Ci P Cj + \/ -Ppi Ppj ff ) sin tpij + Nq 
( V P Ci P Ci + V P D t Pd,H) cos (f>ii + Nj 


( 31 ) 


where N, and Nq are the real and imaginary parts of TV, respectively. Note that although TV/ and 
Nq have zero mean, their joint statistics are still influenced by the correlation between hUt) and n (t) 
These^st atist ics are analyzed in Appendix A, and the density function for the phase estimation error 
A<j>u = <f>u — <fi u is derived. 

In [2], a quantity known as the correlator SNR is introduced, defined as 


SNR corr 


E[Z]E*[Z] 
E\ZZ*} -E[Z]E*[Z] 


(32) 


The correlator SNR is 
to the variance of the 


a measure of the spread of the phase error density p^,(A0 H ) and is inversely related 
phase error. In [1], where FSC is analyzed for independent noises, it is shown that 
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the phase error density can be expressed solely in terms of the correlator SNR For the correlated noise 
case the density is given in Appendix A in terms of the correlator SNR and the correlation parameters 

pu and ipU' 

Figures 7 through 9 show the density function p*(A <f>) for various values of p and f The signal 
parameters chosen for these curves are ( Pt/N 0 )i = (Pr/Noh = 25 dB-Hz, A - 9° eg, wit 
subcarrier harmonics included in the correlation. The correlator parameters are Bj p = B, p - 15 kHz 
and T =3 s. Note that even for a noise correlation as high as 0.4, the density function looks remarkably 
like thatof the uncorrelated noise case. Simulations were performed for the same parameters and densities 
collected for the measured phase estimates. These results are shown with the analytical curves in Fig. . 



4> (deg) 


Fig. 7. Phase estimate density. 

2. Arrayed Symbol SNR and Symbol SNR Degradation. Using the set of estimated phases to 
align the signals, the combined signal becomes 


Vcombit) = Scomb(t) + n co mb(,t) 

L l 

= e i * u s i (t) + 'Y,0i e }4,u n t {t) 


= ypi (y/Pci-iVP^ sqr M + <u) + E 

The combined signal power conditioned on the set of phase errors A<p u is thus given by 

P' T = E [SeombW] E [s*omb(t)} 


(33) 

(34) 

(35) 


L L 




1=1 i=l j-i 


(36) 

(37) 


Similarly, the conditional noise power spectral density is given by 
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Fig. 8. Phase estimate densities: (a) p = 0.4, *F = 0 deg, (b) p = 0 4 4* = 45 dea (cl n - n i 
V = 90 deg, (d) p = 0.4, V. 135 deg, and (I) p = 0.4,% = 1M deg 9 ’ (C) P ~ °- 4 ’ 


2 B ^ ^ lcomb ^)™com.b{t)\ 

(ifil jY^Pij e^ iJ 

Taking the ratio of Eq. (37) to Eq. (39) yields the conditional P T /N 0 of the combined signal, i.e. 

(Pr\' = Pt i gji 7, 2 + gti Ef;; y^e^.-A^ ) 

\N 0 ) N 0 i £\ =1 -y t + ilil ]) 1 ? 2 Pij ei^'j 

After carrier and subcarrier demodulation and matched filtering, the conditional symbol SNR 
arrayed signal is given by 




i= 1 i=i 


(38) 

(39) 

(40) 
of the 
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Fig. 9. Phase «^tedensitle8: o (a) p ;°|^° a ^ e { b p > £ o = .8°y l 1 = M «£?’ ^ P = °' 8 ’ 


SNR' = 


2Pj\ 


Ef., 1? + Ti., Ea: **• 




JV, 


ol 


^ym 7i 


+ Zti E*-> (7i7j) 1/2 P‘j eJ ^’ J e^--^ 


cldc 


c sc^sy 


(41) 


where C c ,C sc , and C ay are the carrier, subcarrier, and symbol reduction functions, respectively. The 
unconditional symbol SNR issued ph^ 

^ — ;" 1 to each of the * 

quantities separately yields an expression for the unconditional symbol SNR, namely 


SNR 


2Pri ci ci d 


N 0 iR. 


sym 


7T 7T 

/-/ 


Eti7? + EtiEfew- 


,j(A0i« - A^jj) 


Ef=! 7. + Ef-1 E# (7i7,) 1/2 P« ^ 


p(A<t>u) ■ ' -p{A<Pil) 


x dA0i2 • • ■ dA<t>iL 


( 42 ) 
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pM 



pM 



PM 


pM 



4> (Peg) 


Fig. 10. Phase estimate densities with simulation points: (a) p = 0.0, (b)p =08 «/ -Odea 
(c) p = 0.8, y = 45 deg, and (d) p = 0.8, V = 90 deg ’ V 89, 

tt e Meli e SNR Si F ^° nS "T SiV6n in Appendix A ' Fin ally. taking the ratio of Eq. (42) to 

the ideal SNR, Eq. (25), yields the degradation for full-spectrum combining: 

D fsc =clcT c c% 


7T 7 r 

“/■■■/ 


7t 2 + 1 

Etl 7 . + Ef=! ^ ei(^^) P{A<i>u} ' ' -P (A ^) 


x dA<f>i2- -dA<t>i L G~. 1 


(43) 


Note that D /sc is equal to one in the upper limit, where A 0 H = 0 for i = 2 ■ ■ ■ L and ( 7 ^ = < 72 " = C T = 1 

sSe ”; 1 , ***■ md 55 a »" f «■» '-P 

C. Simulation Results 

A simple two-antenna array was simulated under conditions of correlated noise to verify the analysis 
given above. The symbol SNR of the combined data was measured using the split-symbol 

Sd m «,e Pn /N = P - S ,raR 0 ' H SNR B t0 °'*‘ in mea ‘ i ' ,red The signal parameters 

rpi . Tl \ 01 . t 2 /No 2 25 dB-Hz, R sym - 200 symbols per second (sps), and A = 90 deg 
he carrier, subcarrier, and symbol loops were operated with bandwidths of 3.5, 0 75 and 0 15 Hz 
respectively, with a symbol window of 1/2. The correlation coefficient between the noises, p, and the 
relative noise phase, were varied over a range of values. 

tor SNR e Th^t° W Vi. mUl ? i0 !i a '°! lg With CUrV6S d6SCribing anal ytical results for a “high” correla- 

SNR. The correlation bandwidths and integration time were chosen so that degradation resulting from 
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Fig. 11. FSC degradation, high correlator 
SNR— theory and simulation. 

imperfect phasing is negligible compared to the carrier, subcarrier, “ d 

that more degradation is incurred with increasing noise correlation for ^ - 0 deg, and that degraaat 
decreases as p increases for = 180 deg. 4 This can be explained by noting the effert ° fv * ry ^ P ^ in 
the arraying gain. For V> = 0 deg, increasing p causes a decrease in arrayed symbol SNR, as explained 
Section III. A. The loop SNR of the three loops, therefore, decreases, resulting in more carrier, > 

A u 1 ™ R V contrast when xb = 180 deg, increasing p increases the combined Pt/N q and raises 

the thnl loop SNrL This results in less degradation in demodulating the signal. Since the 

SNR «gh in this example, the demodulation losses are the dominant source of degradation, and the 

trend shown in Fig. 11 is thus explained. 

Figure 12 shows the same results performed for a relatively “low” correlator SNR. Here, the degradation 
curve for J = 180 deg actually lies below the curve for $ = 0 deg. This result, although seemingly counter- 
intuitive tan nevertheless be explained qualitatively. Note from Eq. (41) that the phase error terms A<fiu 

decreSe the arrayed signal power, but can decrease or increase the arrayed 1 norse power, dqxrtug .on 

* he resuhs ,1 6 m« LmlLTleliation. Therefore, estimating the phase 

^"uL “ a SoSd y penai,y: The combined signal power is les^ned, and the combined „0« 
power increases. This results in increased degradation due to imperfect phase alignme . _ 
hand when xb = 0 deg, phase misalignment decreases the arrayed noise power Since 0 U - 01, “ 
cZ aTiglg the s^s imperfectly also lessens the constructive addition of the noise. The reduced 
noise power due to phasing errors, therefore, has a mitigating effect on the degradation incurred. 

It should be noted that the fact that the y = 180-deg case has more degradation than *** * = 0 ' de ® 
case in this example does not mean that the overall performance of the array is worse for lA 180 «’ 

Recall that degradation is defined as the deviation from the ideal arraying gain, Ga ■ In the above exa p , 


^Tphrase “decreasing degradation” is used loosely to mean decreased synchronization losses; in actuality, numerically 
lower degradation implies greater losses incurred. 
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— a-— -1.VII, IWfT VUII 

SNR— theory and simulation. 

?f* ly h ? her ' the ; dMl sai " ls ■»■— <■* ■*- *-“» “ 
snr, l., ju r “™ y * — * — — 


IV. Complex-Symbol Combining Performance 


A block diagram of the complex-symbol combining arraying scheme is shown in Fi<r n p u • , 

complex-symbol combining are ‘ °T “T ^ 

r ;r s r^ p i~ sr “ d 

Xe a S A^bZS^t; ^'finally ^t’d^od^.tct^ ” “* 

^5H5SS“t/^{£ r co™p«« f or^b“^ 

ZF&z SESSSHS 

A. Antenna Phasing 

The complex-symbol stream from the ith antenna is given by 


— V -^DiCsc iC S y t d(k) ^^!>T a k+d,) 


( 44 ) 
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Fig. 13 . Complex-symbol combining. 

, r C are the subcarrier and symbol reduction functions for the ith receiver, T s is the 

r^boUtae! Ilf*) is the noise output from the ith mulched filter. Taking the complex product 
between the 1st and ith streams yields 


Z - Y\{k)Y*{k) = \/ PdiPd.Csci C aCi C syx C syi + N s 


■Ni(k)N:(k) 


. • 1 • .orm AT (k) has zero mean. Once again, the complex-noise product N\(k)N*{k) 

where the signal-noise term N a , n { ) introduces a bias to the signal correlation 

has nonzero mean if the correlation coefficient is nonzero d g oes not contain 

vector. Note, however, that the spectrum o e signa ^ subcarrier collapses all the data 

empty bands as in the case of full-spectrum combining, ^"“^idth Sin“ t“e shaded information 
sidebands to baseband, allowing a muc narrower com i ^ d t measure the 

rate for CSC is equal to the symbol rate there is no ««ss ^ ^ fo[ ^ 

SCS this possibility, however, we calculate the expectation 

of the noise product, E\Ni(k)N"{k)}. 

r Her the block diagram of Fig. 14, which shows the processing for complex-symbol combining up 
th^matched filter ot^ The signal s,(.) is the subcarrier reference from the rth subcarrier ioop, 

given by 

s t (t) = sqr(o> 3C t 4- 0 9C 4- 4>sa) 

where 6 is the instantaneous subcarrier phase and is the instantaneous phase error in the ith loop, 
for i = l! • , The limits of integration for the ith matched filter are given y 

t u = kT s + Ti (47) 


tu, — (k + 1 )Ts + Ti 


228 














N L (k) 


Fig. 14. Matched filter noise outputs for CSC. 

where Ti is the timing error in the ith symbol loop. The matched-filter 
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noise samples are, therefore, given 
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( 49 ) 


( 50 ) 


be Sr s “ “ d symW *“"* «■» - ti*. 

baseband noises ie® ^P^smns with the cross-correlation function for the complex 

R h„h, ( u,v ) = E [ni(u)n*(v)] = a ly 6(u - v) ( 51 ) 

yielding 
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where the limits of integration of v are given by 

tmin = ma x{kT, + r t , kT s + Tj) 

t m ax = min {{k + 1 )T» + t u ( k + l)T, + Tj) 


(53) 

(54) 


Finally, integrating with respect to v yields 
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Note that, i„ the absence of phase errors in an, of the kjW B* oJ 

- «*«■ *• 

Sctlr symbol phase errors also reduce the noise correlstron rt tins pent. 

Calculating the unconditional romance of 

of Eq. (57) with respect to the phase errors .0.**. assumed to be Gaussian. This 

perform this computation. F ' rst ’ 1 6 10 dB and is consistent with the approximation made in 

condition is nearly satisfied for loop mutually independent. This statement is not 

[1], Second, the phase errors of all loops are are affected by the same 

strictly justifiable, since the subcarner an sy Darate receivers are correlated. Nevertheless, it is 

noise and, furthermore, because the noises view y P evaluating the unconditional covariance. 

invoked for the purpose of making a ^^^ ^Gaussian-distributed with known mean and variance, 
The quantities <j> SCi - <t>sc, and <p sy , - <Psy, are uien « 

and the unconditional expectation E[Ni{k)N*{k)} becomes 


E[Ni(k)N*{k)] 


= aae 
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(58) 


Equations (58) and (45) can be combined to calculate the ratio of the signal-t^noise correlation magni- 
tude, analogous to that computed in (29): 
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noise correlation is a significant bias in estimating the relative signal phase. If |5|/|N| is much less than 
I, then an extra correlation is needed to compensate for the noise vector, as mentioned earlier. On the 
other hand, if this quantity is much greater than 1, then it is unnecessary to add the extra matched filter 
channel to perform the noise-only correlation. Note that collapsing all the data sidebands to baseband 
and performing matched filtering before the correlation takes place substantially decreases the correlation 
bandwidth relative to that of the FSC scheme described in [2], The full-spectrum combining scheme can 
optionally be modified to employ a similar strategy by using a series of matched filters for each sub carrier 
harmonic, as discussed earlier. Estimating the degree of correlation p that will be observed for a particular 
antenna pair and applying the rule described above will indicate whether or not the noise contribution to 
the total correlation is substantial and must be compensated for by performing an additional correlation. 

Here we briefly describe how the extra matched-filter outputs can be used to measure the noise cor- 
relation: The complex baseband signal from each antenna can be shifted in frequency so that an empty 
portion of the spectrum is located at baseband. This may be accomplished by shifting by an even multiple 
of the subcarrier frequency, i.e., 

m) = +jy/P^d(t) sqr(w sc t + 0 ac )e^ t+ ^ + #*(«)) e? Nu -'* = s'^t) 4- n'(t) (60) 

where N is an even integer. The shifted signal can then be multiplied by the subcarrier reference from 
the ith antenna and passed through a matched filter using timing from the zth symbol loop, as shown in 
Fig. 14. Thus, 


1 r(k+l)7 1 s+ r i 

N'(t) = — / nj(i) sqr(a; 5C * + 6 SC + <p sci ) dt (61) 

JkT s +n 

From the above analysis, it is clear that E[Nl(k)Nj(k)] will be given by Eq. (58). Correlating the two 
noise-only matched filter outputs then yields a quantity that can be subtracted from the total correlation, 
Z, to compensate for the noise bias. The density function for the phase estimate computed using this 
technique is similar to the FSC case and is analyzed in Appendix B. Note, however, that performing 
this compensation requires increasing the combining bandwidth beyond what is required for CSC in the 
uncorrelated noise case, as well as additional hardware to process the extra channel containing noise only. 
A tradeoff in performance versus complexity must, therefore, be made to determine if comp lex- symbol 
combining is an attractive option when correlated noise is present. 

B. Arrayed Symbol SNR and Symbol SNR Degradation 

An expression for the conditional arrayed symbol SNR can be obtained in a similar manner as is the 
full-spectrum combining case. The combined signal for complex-symbol combining is given by 

Y comb (k) = S comb (k ) + N comb (k) = (V^Csc i C sy Ak)e iiu,bT ’ k+0i) + Ni(k)) (62) 


The conditional signal power, defined as E[S com b(k)]E[S* omb (k)}, is given by 


P zomb — Pd x 
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(63) 
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where, as before, A <pu is defined as the error in estimating the phase difference between the 1st and ith 
signal, 4>u - <f>u. The one-sided power spectral density of the real and imaginary parts of iV com &(fc) is 
given by 


N 0 = T s Var (w com6 (fc)) - T s E 


1 / \J = 1 


(64) 


Using the relations 
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Pii \/No No e^> 

EmvN-ik)} = PiiV ° T °’ — 


Csc,, Csyij 


(65) 

( 66 ) 


Eq. (64) can be shown to be equal to 


K = N 01 


L L 




i=l i= 1 J’ = 1 


(67) 


Taking the ratio of Eq. (63) to Eq. (67) then yields the combined P D /N 0 for CSC. The combined signal 
is finally processed by a baseband Costas loop, and the conditional SNR adding in carrier losses is given 

by 


SNR' 


2.P j 


D, 


1 (J 


N 0 , Rsym 


Ef=i 7i + Ef=i £f-« 


(68) 


Computing the unconditional symbol SNR requires taking the expectation of the above quantity with 
respect to the phase errors (f> SCi and <f> 9y . for i = 1, • * • , L, the phase estimates <j>n for i = 2, • , L, and 

the carrier phase error <f) c . Once again, we assume all loop phase errors and phase-aligning errors are 
mutually independent. Thus, integration over the carrier phase error <f> c is accomplished easily by consid- 
ering the carrier reduction function C\ separately. However, unlike the case of full-spectrum combining, 
the subcarrier and symbol phase errors appear in both the numerator and the denominator. The expec- 
tation with respect to the subcarrier and symbol phase errors, therefore, cannot be given in closed form. 
Calculating the unconditional symbol SNR for even a simple two-element array would thus require a fifth- 
order numerical integration. Rather than resort to such brute-force tactics, we make further simplifying 
assumptions to allow evaluation of some of the integrals in closed form. 

In taking the expectation with respect to the 4> 9Ci and <j) syi terms, we apply the approximation 
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to the ratio of Eq. (68), yielding 
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where * 4C is the set of subcarrier phase errors cj> aCx for i = 1, • • • , L, is the set of symbol phase errors 

Fa i<f 1 ’ ’" w ’ $ 18 the ; S6t ° f phase estima tes for i = 2, • • ■ , L. The approximation of 

Eq. (69) is reasonable if the mean of y squared is much greater than the variance of y (i.e., if y is nearly 

1S met /° r ^ , C r e Under COnSideration ' since it is implicitly assumed that 
tvni u mu SUb< j arner and symbo1 lo °PS are high enough to maintain lock, with 13 dB being a 

vZe e ' h ’ n VarianC f S 0f ? e redUCti ° n fUnCtl ° nS and C *»» which contain the loop 

phase errors, will be small compared to the mean of the entire denominator term. 

By the above argument, the unconditional SNR can be evaluated as 


SNR = 
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(71) 


since^N/? ?, n ^ for complex-symbol combining is identical to that for full-spectrum combining; 
sine e SNR ]d t is defined as the SNR that would be obtained in the absence of synchronization errors 
is value is independent of the order in which combining and demodulation occur. Thus, the degradation 
complex-symbol combining is found by combining the results of Eq. (71) with Eq. (25), yielding 
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(72) 


C. Simulation Results 

Simulations of a two-antenna complex-symbol combining system were performed. The signal pa- 
rameters^used were the same as those used for the full-spectrum combining simulations: P T JN 0 = 
tJ 02 5 dB-Hz, R syrn - 200 sps, and A = 90 deg. The loop bandwidths were also set as before; the 
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carrier subcarrier, and symbol loop bandwidths were 3.5, 0.75, and 0.15 Hz, respectively with a symbol 
window of 1/2. Both the compensating and noncompensating methods of estimating the signa P 
difference were implemented. In Figs. 15 and 16, simulated and analytical degradation values are 

for various values of p and ip. 

For the uncompensated case, the degradation curve drops down sharply for 0 = 90 deg and 0 = 
180 deg. One cause for this is the bias in the complex correlation used to estimate the relative sign 
phase For the parameters being used, \S\/\N\, given by Eq. (59), is equal to 3.15 for p = =0.5. Thus, the 
Lise vector is of comparable but lesser magnitude to that of the signal in .the Note 

that for 0 = 0 deg, the noise correlation phase is equal to the relative signal p (0 0), 

vectors S and N are colinear (see Fig. 5). The noise vector, therefore, does not bias the measuremen 
away from the desired quantity, and the downward trend is not present. 


For the compensated case, less overall degradation is observed. However the 0 = 180-deg curve still 
droos down with increasing p. Recall from Section IV.A that imperfect subcarrier and symbol tracking 
tend to decrease the power levels of the individual signals at the matched filter output and decrease the 
matched filter noises. When 0 = 0 deg, this has a beneficial eff^ ^n the ^ 
SNR since it reduces the coherent addition of the noise. By contrast, when 0 - 180 deg a high degree 
correlation between the noises is desirable, so that the noise cancels maximally. Thus decreasing 
this correlation lessens the arrayed SNR and causes more degradation. This explains the fact that 



Fig 15. CSC degradation, phase uncompensated: 

theory and simulation. 


Fig. 16. CSC degradation, phase compensated: 
theory and simulation. 


V. Example: Galileo Scenario 


I„ order to illustrate the major concepts presented in this article, the perform^ of Ml-spectrum 

combining and ^ ““^iency '(HEF) 'amenna is ciin"n for this example. 

ZXt LaS strength and degree of noi* correlation are 
developed- These quantities are then used to calculate the arraying gam and degradation for each of 
two schemes. 
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A. Signal Parameters 

Of h\ f. he C f 3Se ° f th A 6 G , alile ° spacecraft ’ correlated noise will be contributed by Jupiter being in the beam 
of both antennas. As discussed in Section II, the contribution of a background body to total system noise 
depends on its angular separation from the spacecraft and on its total flux, which varies with its distance 
rom Earth. Values for the Jupiter-Earth probe (JEP) angle and Jupite^Earth 
m ephemens information for the Galileo tour. For the purpose of this example, we select values that 
maximize the noise contribution of the planet to estimate the impact of correlated noise in a worst-case 
scenario. Thus, we assume the JEP angle is zero and that the Jupiter-Earth range is at its minLum 

fo/oSS^and DSS 15 * 'r^ ~ fin U f mg theS6 ValueS ’ the temperature contribution of Jupiter 
for DSS H and DSS 15 are T s , = 6.6 K and T, 2 = 1.4 K, respectively. Note that the temperature 

contribution is higher for DSS 14 due to the greater aperture size and antenna efficiency. 

fnrTh 6 7 n ediC n e lf gnal parameters are 35 follows: (Pt/N 0 ) i = 22.0 dB-Hz and {P T /N 0 ) 2 = 11.6 dB-Hz 
for the 70- and 34-m antennas, respectively; A = 90 deg- and R - 200 sns Nnfc rl, e ■ 
awiimin^ iho ^ u n *vm - sps. Note that since we are 

strenrth 12 ^ s P a cecraft are at their closest range, the spacecraft signal is also at its peak 

DSs S ^nd DSS lT ° f The t0tal SyStem tempera tures predicted for 

iJoo 14 and Dbo 15 are 22.6 and 42.2 K, respectively. 5 

To determine the degree to which the source is resolved on this array baseline, we must compare the 
ringe spacing to the angular size of the source. In our example, the observing frequency f 0 is 2.3 x 10 9 Hz 
and the maximum possible projected baseline is the physical separation between the two antennas which 
s approximately 500 m. Thus, the smallest possible fringe spacing is 2.5 x 10~ 4 rad. At a range of 4 0 AU 
Jupiter has an angular size on the order of 1 x HT 3 rad. Since these values are comparable we cannot 
use either the long baseline limit or the short baseline limit in evaluating p (see Section I). However for 
e purpose of determining the impact of the correlated noise in the most extreme case, we overestimate 
the degree of noise correlation using the upper bound on p, given by 


P = 


T tl T a , 

TyT 2 


■ 0.1 


(73) 


B. Arraying Performance 

Using the two Pj/N 0 levels and correlation coefficient p found above, the ideal arraying gain G A can 

Not Tfr, " 3 UnCtl ° n ° f * Eq - (25) ‘ A gmph sh0win S this relationship is shown in Fig 17 

Note that the arraying gain in this example is much smaller compared to our previous examples of two 

q a antennas, since the signal level of one antenna is approximately 10 dB lower than the other. For 

Bv g ;Imr^t re h t T P ° ne f ° f the noise adds maximally in phase, thus decreasing the arraying 
gam. By contrast, the background noise interferes destructively for V- = 180 deg, resulting in greater 

rsr- r e ? correlation coefficient is re,ativeiy iow in this examp,e > 

tne best-case and worst-case scenarios is only about 0.45 dB. 

anH R 0 P 07 e H ntatiVe ^ the f “ nfer . subcarrier, and symbol loop bandwidths were chosen as 1.5, 0.4, 

and 0.07 Hz, respectively. For full-spectrum combining, a correlation bandwidth of B COTr = 2 kHz wai 
used, with a correlation time of 15 s. The total degradation for FSC as a function of * is shown in Fig 18 
ong with simulation points. Because the correlation coefficient p is relatively low in this exampk the 
degradation is almost constant with respect to the phase parameter The combined P T /N 0 only varies 

chan° ff e S m y h d ? “ t rangeS . from 0 to 180 de & thus > the loop SNRs of the three loops also do not 
change much, and synchronization losses remain essentially constant. 


5 Predicts for noise and signal parameters were obtained from the Galileo S-Band Analysis Program (GSAP). 
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Fig. 18. FSC degradation for Galileo 
signal parameters. 


The same signal parameters and loop bandwidths were used to simulate the 
case A slight variation of the basic scheme, known as complex-symbol combining .with aiding (CS ), 
was implemented. This scheme is discussed in [2] as an option for arraying the Galileo signal In i CSCA, 
the subcarrier and symbol references from the receiver tracking the stronger signal are used to track 
signal from the 34-m antenna as well. This technique can be used to perform complex-symbol com mi g 
even if the 34-m antenna signal is too weak to achieve subcarrier and symbol lock on its own. Thus, 
the loop SNRs for the 34-m antenna subcarrier and symbol loops are equal to the correspon g 

antenna loop SNRs. 


Equation (59) can be applied to determine whether or not the “noise-only” channel is needed to phase 
the array. Substituting in values from above, we find 
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= 2.39 


(75) 


Thus the magnitude of the noise correlation vector is less than but comparable to that of the signal 
correlation vector To illustrate the impact of the phase bias in aligning the signals, CSCA was simulated 
wTh botl^ ^compensating and uncompensating method for estimating the relative signal phase. In 
Fif 19 we showThe degradation for CSCA for these two cases. The correlation time used to estimate 

the relative signal phase was 2 s. Note that a shorter estimation interval than the full-spectrumcombming 
the relative s fT p effectlve corre i ation bandwidth is equal to the data bandwidth of 200 Hz 

aTop C posed to 2 kHz for FSC. For the compensated case, the degradation is essentially constant since, 
as opposed to 2 does not affec t synchronization losses much. For the uncompensated 

r the degradation becomes greater as the difference between the noise and signal phase t f> grows larger, 
since the noise correlation begins to bias the phase estimate further away from the relative signal phase^ 
?hL effect can be seen graphically by referring once again to Fig. 5, where the complex-signal and noise 

correlations are represented as vectors. 
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Fig. 19. CSC degradation for Galileo signal parameters: (a) phase compensated and 

(b) phase uncompensated. 


VI. Conclusion 

The effects of correlated noise on the full-spectrum combining and complex-symbol combining arraying 
schemes have been analyzed. As seen in Section II, accurate modeling of the noise correlation properties 
for a given antenna pair requires detailed analysis of factors such as the source structure and position, 
the antenna gain patterns, and the geometry of the array. However, the correlation coefficient can be 
determined easily in cases where the baseline is either very short or very long. These two extreme cases 
can be used to obtain a rough idea of what degree of noise correlation can be expected for a given scenario. 

Describing the correlation between the various antenna pairs in an array by the parameters pij and 
expressions for the ideal arraying gain and arraying degradation were derived. Several important 
differences from the uncorrelated noise case were noted. For a given set of signal levels (. P T JN 0i ), the 
ideal arraying gain when the noise is correlated may be higher or lower than when the noise waveforms 
are independent. This reflects the fact that the noise may add constructively or destructively, depending 
on the relative signal and noise phases (i.e., the ipij parameters). 

In addition, correlated noise can have a significant impact on the synchronization processes used to 
combine and demodulate the signals, which vary with the specific arraying technique used. Most notably, 
a bias due to the noise correlation is present in the conventional method of estimating the relative signal 
phases. Since the magnitude of this bias is proportional to the correlation bandwidth used, full-spectrum 
combining is potentially more sensitive to this problem than complex-symbol combining, depending on 
the specific method used to correlate the signals. A modified method of phase estimation, where the 
correlation due to the noise alone is measured and compensated for, can optionally be employed for both 
FSC and CSC, a s necessary. 
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Appendix A 

Performance of the FSC Correlator 


For full-spectrum combining, the phase difference between two signals is estimated by performing one 
lowpass and one bandpass correlation, as described in Section III.B. After being filtered to some lowpass 
bandwidth, Bi v Hz, the signals from antenna 1 and antenna i are given by 
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where the subcarrier is expressed in terms of its sinusoidal components that are passed by the lowpass 
filter. The two signals passed through the bandpass filter of bandpass B(, p Hz contain only noise and are 
given by 


Vbpi W “ ^6pi (0 


(A-3) 
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The complex quantity used to estimate the relative signal phase <p u = 6 X - is given by 


Z = I + jQ 


2 r B 1 f 

~~ y* / Vipi ^ ~ T> 7^ / ^pi Vbp t dt 

x corr J £>bp J- corr J 


™bpi ^ . dt 


~ (VP Cl Pc, + V-PdjPd.h) e^ 1 - + f (n s n + h lpi n^ p% ) dt - ^ -L- f j 

corr J P^bp corr J 

= (VP^Pci + y/P&ZH) e?*" + TV (A . 5) 


In most cases, the contribution of the signal-noise term n s , n (t) to the total noise power is much smaller 
than that of the nois^noise terms, and can be ignored. This is especially true if the P T /N 0 levels of the 
two signals are very low, or if large correlation bandwidths are used. By the Central Limit Theorem, 
the complex noise N can be approximated as Gaussian if the correlation extends over many independent 
samples (i.e., if T corr is much greater than the inverse correlation bandwidths). After averaging, the 
variance of the real and imaginary parts of N can be shown to be equal to 


A/ = Var(Nj) = ( B lp + ^ ) (N 0l N 0i + a 2 u cos20£) 

1 corr \ &bp ) 


(A-6) 


1 /- K 


\q = Var(N Q ) = B ip + (N 0l N 0t - a 2 u cos2^ t ) 

1 corr V t>bp J 


(A-7) 


where N, and N Q are the real and imaginary parts of N, respectively. The covariance of N { and No can 
be shown to be equal to v 


A/q - Cov(Ni, Nq) - — \Bi p + sin 2^ (A-8) 

Furthermore, it is clear from Eq. (A-5) that the means of the real and imaginary parts of Z are given by 

m i ~ ( yj ■ P Cl P Ci + \JPd x Pd x cos <f> u (A-9) 

m Q = {y/PcxPci + y/PoxPDiH^ sin i (A-10) 

Equations (A-6), (A-7), (A-9), and (A-10) can be combined to compute the correlator SNR as defined in 
[1], i.e., 
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SNR 


corr,fac 


E[Z]E*[Z] 
E[ZZ*\ -E\Z)E*[Z) 


mj + 

A / -f A q 


2{B lp + {B? p /B bp )) 



(A-ll) 


Equations (A-9), (A- 10), and (A-6) through (A-8) can be used to determine the joint density function 
p IiQ (I,Q). Since the density of 4> li = tan -\Q/I) is the desired quantity, we express the joint density 
function in terms of polar coordinates, using the variable definitions 


r A ^72 + Q2 


(A-12) 


0 = tan" 



(A-13) 


The density function for jointly Gaussian random variables is given in polar form by 


/ r ,^( G < A ) 2tt(\i\q - A iq) 


x exp 


\i(r cos <j> — mi) 2 


2A/q(rcos<ft — m/)(rsin^> — mo) + Aq(rsin<ft mo) 
2(A/A q - Xiq) 


(A-14) 


Integrating Eq. (A-14) with respect to r yields the marginal density of </> alone. Expressing the phase 
estimate density in terms of the estimation error No = cpn - <pu yields 


U(Aj>) = Cl exp (^SNRcorrJsc - ^4 —) l 1 + S* G * e ° l erf Ga + 1} ] (A ' 15) 


where 


1 - p 4 

Gl 27 t( 1 - p 2 cos(2V> - A0)) 


(A-16) 


/rii> — cos A0 - p 2 cos(2V> - A<£) 

G 2 = y'SNRcorrjsc (1 _ p4 )(! _ p2 cos (2i/> _ A0)) 


(A-17) 
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Appendix B 

Performance of the CSC Correlator 


The method of estimating the relative signal phases for complex-symbol combining is analogous to 
the full-spectrum combining algorithm; using the extra correlation to compensate for the noise bias, the 
complex correlation can be expressed as 


k=i v *:=i 

i n 

= y/P Dl Pa, c aci C ayi C SCi + — £ y/T%C aCl C ayi ei°'N?(k) 


+ E '/P^C ac ,C ay ,e~^N' l (k) + jj'£N l (k)N;(k) - j- 

k= 1 k = 1 V fc=l 

= >/ -Pdi Pd x C SCl C syi C SCt C ay .e?* u + N (B-l) 

where N is the number of symbols averaged over, given by TV = T corr /T sym , and the noise term N has 
zero mean. The statistics of this noise can be analyzed in the same manner as before; here, the effective 
correlation bandwidth for both the lowpass and the bandpass correlation is R sym /2. Using the definition 
given by Eq. (32), the correlator SNR can be shown to be equal to 


SNR 


Pd i 


T C orrC sc\ 


■corr.csc 


i 


'*Vi 


N n 


Cl e, Cj v% + Cj ci C? yi ( l/7i) + {NJP Di )2R 


(B-2) 


aym 


The density function for the phase estimation error can be found in a manner analogous to that applied 
in Appendix A. The only difference is in the expression for the correlator SNR; otherwise, both problems 
are inherently governed by the same mathematics. The density function for the phase estimation error 
A <j>u is thus given by Eq. (A-15), with SNR corr j sc replaced by SNR corryCSC . 
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This article presents a study of seismic data compression techniques and a com- 
pression algorithm based on subband coding . The algorithm includes three stages: 
a decorrelation stage , a quantization stage that introduces a controlled amount of 
distortion to allow for high compression ratios, and a lossless entropy coding stage 
based on a simple but efficient arithmetic coding method. Subband coding methods 
are particularly suited to the decorrelation of nonstationary processes such as seis- 
mic events. Adaptivity to the nonstationary behavior of the waveform is achieved 
by dividing the data into separate blocks that are encoded separately with an adap- 
tive arithmetic encoder. This is done with high efficiency due to the low overhead 
introduced by the arithmetic encoder in specifying its parameters. The technique 
could be used as a progressive transmission system , where successive refinements of 
the data can be requested by the user. This allows seismologists to first examine 
a coarse version of waveforms with minima] usage of the channel and then decide 
where refinements are required. Rate-distortion performance results are presented 
and comparisons are made with two block transform methods. 


I. Introduction 

A typical seismic analysis scenario involves collection of data by an array of seismometers, transmission 
over a channel offering limited data rate, and storage of data for analysis. Seismic data analysis is 
performed for monitoring earthquakes and for planetary exploration, as in the planned study of seismic 
events on Mars. Seismic data compression systems are required to cope with the transmission of vast 
amounts of data over constrained channels and must be able to accurately reproduce both low-energy 
seismic signals and occasional high-energy seismic events. 

We describe a compression algorithm that includes three stages: a decorrelation stage based on subband 
coding, a uniform quantization stage, and a lossless entropy coding stage based on arithmetic coding. 
Rate-distortion performance results are presented and comparisons are made with two block transform 
methods: the discrete cosine transform (DCT) and the Walsh-Hadamard transform (WHT). 

Subband coding methods are particularly suited to the decorrelation of nonstationary processes such as 
seismic events. For most seismic data, signal energy is more concentrated in the low-frequency subbands, 
which suggests the use of nonuniform subband decomposition. The decorrelation stage is implemented 
by quadrature mirror filters using a lattice structure. Adaptivity to the nonstationary behavior of the 
waveform is achieved by dividing the data into blocks that are separately encoded. 
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